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Abstract 

A mathematical model of the distribution function for the discrete 3- 
disk is proposed in order to utilize in the statistical evolution equation of 
the 3-dimensional Universe. The model is constructed based on analyses 
in known exact solutions of recursion equations for the generating func¬ 
tions of the discrete 2-disk. The model distribution exhibits three types 
of phases characterized by geometrical structure of the disk with either 1, 
2, or 3-dimensional nature. Transitions between those phases are either 
cross-over, 1st order, or 2nd order depending oil the model parameters. 
The proposed distribution function is compared with numerical simula¬ 
tions of the dynamical triangulation(DT) with S 3 , and D 3 topologies. 

PACS numbers: 98.80.Qc 02.50.Ga, 04.60.-m, 98.80.Bp 


1 Introduction 

Discretization of space is one of the standard techniques for quantizing fields 
non-perturbatively. In case of the gravitational field a computational method, 
known as the dynamical triangulation(DT)|I i 21 [3, has been popularly tried 
up to four dimensions. In 2-dimension numerical results have satisfactorily re¬ 
produced analytical expectations. For an example, the string susceptibility pre¬ 
dicted by the Liouvillc field theory[3], [5] E] was obtained for S 2 and T 2 topologies 
with several central charges within statistical errors [7]. Great many numerical 
efforts have also been poured on higher dimensional cases |Hj. Those investiga¬ 
tions have given some insights on the possibility of quantizing the gravitational 
field in four dimensions. For an example, the grand-canonical DT simulation{9] 
has shown the existence of second order phase transition when space is coupled 
to a certain number of matter fields. However, no conclusive results were drawn 
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at that time. It may be not only because of lack of the computational ability 
but also lack of analytical supports for the method in contrast to the lattice 
QCD. 

While the discretization of space has been originally introduced as a tech¬ 
nical tool to avoid the ultra-violet problem of quantizing gravity, there exists a 
possibility to regard the discrete space as the fundamental structure of space- 
time[ini HU HT2] . Based on this idea we have constructed the ^-dimensional 
Universe as a complex of d-simplices with the d-disk topology, where the S' d_1 
boundary is assigned as space, while the direction perpendicular to the bound¬ 
ary surface as time. Statistical evolution equation of the Universe has been 
formulated as the Markov process jT3], and applied to the 2-point temperature 
correlation observed by the WMAP[I2j in 4-dinrension. In this model the evo¬ 
lution is regarded as a kind of the diffusion limited aggregation of simplices 
derived by entropic pressure. For obtaining the distribution function of sim- 
plicial complexes we had to employ numerical simulations which have been the 
unique machinery beyond 2-dimension up to 4-dimension at that time. 

As far as the 2-dimensional Universe is concerned mathematical studies of 
discrete geometry have shown significant results, such as giving analytic forms 
of distribution functions for various cases of discretization. These results are 
elaborated mainly due to the graph counting theory originated by Tutte and 
his followers jl5] [T51 ilTi . As an important application of the results the master 
equation for evolution of the 2-dimensional Universe has been solved analyti¬ 
cally with using asymptotic limits of the exact distribution functions for the 
force field [20]. However, in higher dimensions no analytic form of the distribu¬ 
tion function was available, and all the essential properties were gathered only 
through numerical simulations, for examples [2T1 12211221 [221 [2511251 [27] . in the 
3-dimensional caseO] 

Recalling numerical efforts without sound theoretical support can conclude 
only a temporal consolation in the history of lattice quantum gravity, we focus 
the main aim of this paper to propose an analytical model of distribution func¬ 
tion for the discrete manifold. We hope it will be able to provide a theoretical 
support for numerical methods and serve an initial input for the statistical evo¬ 
lution equation of the 3-dimensional Universe. As one of the first trial we would 
like to propose an analytic distribution function for the discrete 3-disk which is 
constructed based on analyses of known exact results in two dimensions. 

This paper is organized as follows. In order to appeal the model to be 
reasonable, several distribution functions in 2-dimension are reviewed in the 
following section (Section 2). In Section 3 a model distribution function for the 
discrete 3-disk is proposed and its asymptotic limit is calculated. The predicted 
distribution function is applied for the S 1 * 3 and D 3 manifolds to compare with 
known numerical results obtained by DT simulations in Section 4. The last 
section (Section 5) is devoted for discussions and comments on the evolution of 
the 3-dimensional universe and extension to the 4-dinrensional manifold. 

1 Recently, the colored tensor model has been developed as an analytic model for higher 

dimensional counter part of the matrix model [28]. We shall discuss relationships to our model 

in the last section. 
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2 Two dimensional distributions 


2.1 analytical models of distribution functions 

Let us consider a triangulated 2-dinrensional plane with the disk topology. The 
triangulation is rooted by choosing a vertex on the boundary as the root ver¬ 
tex, and distinguishing the oriented edge on the boundary incident to the root 
vertex as the root edge. The root edge is pointing under the right-hand screw 
rule directing upwards perpendicular to the plane. Vertices and edges are called 
external if they lie on the boundary, and the remaining 1- and 2-simplices are 
internal. The triangulation with n internal and m + 3 external vertices is said 
to be of type [n,m]. Two triangulations T and T' is isomorphic when a one- 
to-one mapping exists from i-cells(* = 0,1, 2) of T onto i-cells of T', preserving 
incidence including root simplices. Enumeration of the number of isomorphism 
classes of rooted triangulations has been initiated by Tutte[T5]. and subsequently 
refined by Brown [16] and Mullin[T7| and many others for various cases of dis¬ 
cretization. 

Tutte considered a class of triangulations in which no internal edge is incident 
with two external vertices (3-connected). The distribution of triangulations of 
type [n, m] is shown to be 


fn, 0 = 2 


1 ) 

3 ^2 ( m + j + 2) (m — 3j) 

3=0 

An + 3m — j + 1 \ / m + 2 \ / m — 1 

n — j — 1,3n + 3m + 3/ \j + 1, m — j + 2/ \ j. m — j 
An + 1 
n + 1, 3n + 2 


(m ^ 0), 


( 1 ) 


where the suffix T represents Tutte. Here, we have defined the generalized 
binomial coefficient by 

( 2 ) 


n i 

n 2 ,n 3 


Till 

n 2 \n 3 \ 


Allowing two end vertices of an internal edge are external ( 2-connected ) 
Brown obtained the distribution function for type [n, m] to be 


f B = o ( 2m + 3 \( An + 2m + 1 \ 

Jn ’ m \m,m + 2J\n,3n + 2m + 3J' 

where the suffix B stands for Brown. These two classes of triangulations are 
required to satisfy that (a) no two edges have the same end vertices, and (b) 
every edge is incident with two distinct vertices. Mullin considered a class of 
triangulations without the constraint (a) (multiple-edges), and obtain 


/ 2m + 3 \ ^ n+ i / 3 n + 2m + 2 \ 

\m + 1, m + 1 J \n, 2n + 2m + AJ ’ 
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where the suffix M represents Mullin. Recently, Krikun [T5] considered the de¬ 
generate triangulations without the constraint (a) and (b) (loops). The distri¬ 
bution function is shown to be 


( 2m + 6 \ 2 3n- 4 f |n + m+i \ 

\to + 2,?ti + 3 / \n, \n + to + |/ ’ 


(5) 


where the suffix K stands for Krikun. Brown[TS] also considered the distribution 
function for the quadrangular dissection of type [n, m\ with n internal and 2m+4 
external vertices. The distribution function is given to be 


s „„ = 3 f 3 ” + 4 V 3 " + 3m + 2 V (6) 

y ’ \m,2m + 3j \n,2n + 3m + 4:J w 

To our knowledge they exhaust all the known exact distribution functions 
of the 2-disk. Analytical methods for extracting the distribution function from 
the generating function are often highly technical and sometimes intuitive, and 
we recommend readers to refer original papers for the details. 

It will be essential in our later construction of the distribution function for 
the discrete 3-disk to keep in mind that distribution functions of the Mullin, 
Brown, Krikun triangulation, and the Brown quadrangulation have the prod¬ 
uct form of two generalized binomial coefficients. The 2-dimensional mani¬ 
folds considered in these four cases have a common feature that there exists 
a 1-dimensional sub-manifold embedded in the 2-dimensional plane. By ex¬ 
tracting n = 0 components of distribution functions of triangulations, we find 
fom IK = Cm+ i, where C m is the Catalan number, 


C„ 


to , m + 1 / 


and the distribution function of quadrangle dissection gives 


90,m 


/ 3m + 2 \ 
\m 2 to + 3/ 


(7) 


They all correspond to the distribution functions of branched trees made of 
either trivalent vertices for triangulations or tetravalent vertices for quadrangu¬ 
lation. In contrast, the Tutte distribution function does not have such a simple 
product property, since it prohibits 1-dimensional sub-graphs. We consider this 
is the main reason of the product nature of distribution function of maps which 
possesses the hierarchical structure. 


2.2 asymptotic distributions 

The asymptotic limit of a distribution function is achieved by taking the total 
number of 2-simplices large. Writing N t and W for i = 0,1 as total numbers of 
the internal and the external i— simplices, respectively, we have the relationship 
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either N 2 = 2Nq + N\ — 2 for triangulations, or N 2 = No + Ai /2 — 1 for 
quadr angulations. 

Let us consider the asymptotic limit of a distribution function taking the 
Mullin distribution function as an example. Other cases are treated exactly in 
the same manner as this example except the Tutte distribution function, which 
will be considered separately. We introduce the new variables N = N? — 1, and 
q = m/N together with p = n/N . Then, we have the relationship, p = 1(1 — g), 
which implies the ranges of these variables as q = [0,1], and p = [0,1/2]. By 
making use of the Stirling formula the generalized binomial coefficient has the 
approximate asymptotic form, 

a\N + b\ \ 

a 2 N + 6 2 , a 3 N+ b 3 ) 

x 

where we write c = b± — 62 — & 3 , which will represent the critical exponent 
associating to a map whose distribution is expressed by the generalized binomial 
coefficient. Applying the form to the Catalan number, where a± = 2, 02 = a 3 = 
1 , b\ = 62 = 0 and b 3 = 1 , the asymptotic limit is written as 

C N -> N~ 3/2 2 2N 


b i + l/2 


,5 a +i/2 / ,6 s +i/2 V a“ 2 a“ 3 


N 


N c ~ 1/2 ( 8 ) 


jy(ai— a 2 -a 3 )N e -(ai-a 2 -a 3 )N 


within an overall constant, where the exponent —3/2 is characteristic to the 
1-dimensional map. We notice that the most diverging factor TV^ 1-02-123 ^ 
disappears in this case. We shall observe that this diverging factor disappears 
in all the known distribution functions in the following analyses, and thus, the 
distribution is bounded exponentially. 

The asymptotic limit of the Mullin distribution function expressed in terms 
of new variables (N, q) can be written as 

f N (q) N~ 2 q 1/2 (f>(q)$(q) N , (9) 


where 


<P(q) 

$(g) 


(3 + qf/ 2 

(l-g) 1 /2(i + g )9/2’ 

0-1/209/2 (3 + g)l 3+ g)/ 2 

(1 — <jr) C 1 —'J ')/ 2 (1 + q) 1+ 1 


( 10 ) 

( 11 ) 


In the asymptotic limit the distribution is dominated by the term $>(q) N . Since 
the second derivative of the entropy density, which is defined by 


sjv(g) 


-^ lo s[/jv(g)], 


( 12 ) 


is given asymptotically as 


4 

(i-g)(i + g)(3 + g)’ 
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it is positive in the allowed range of q = [0,1], and the entropy density is 
a smooth concave function in this range. The free-energy density defined by 
—S]y{q) + M<7i then, has a peak at go in [0,1) for appropriate values of the 
parameter /i = /x 0 , satisfying 


Mo 


ds N {q) 


dq 


§(q) N is well-approximated by the Gaussian function for large N as 


*(qf ~ <f>(qo) N eM-laoN(q-q 0 ) 2 }, 


(13) 


with 


do — 


d 2 s N (q) 


dq' 2 


We obtain the asymptotic limit of the Mullin distribution function as 
fN,m exp(-^-m 2 + fj,m — XN) 


(14) 


within a constant, where fi = aoqo, and A = aoq^/2 — log[<l>(go)]. As for the 
g-dependence of the function </>(g), which is regular in the range q = [0,1). 
We approximate it by replacing q to its peak value at qo, and absorb in the 
overall constant. The exponent —5/2 of the variable N is characteristic to the 
2 -dimensional map. 

The asymptotic distribution functions of the Brown triangulation, the Brown 
quadrangulation, and the Kulikun triangulation can be obtained precisely in 


the same form as eq.(14) with the same exponent, but different non-universal 
constants (ao,/j,,X). 

As we can expect from the non-product form of the Tutte distribution func¬ 
tion, derivation of its asymptotic limit is a little involved comparing to the 
previous cases. Nonetheless, the method is basically the same as others. We 
replace the intermediate sum over j by the integral over x = j /to, regarding the 
distribution to be a continuous function of x assuming m(= qN) is large. The 
asymptotic form of the Tutte distribution function is then expressed as 


f T N {q)^N~i 


dx (j) T {x,q)<b T {x,q) 


N 


(15) 


within a constant, where 

<h( x >q) = 

and where 

®T{x,q) = 


{2 + (1 — ir)g} 3 / 2 {l — (1 + 2x)q} 1 / 2 
(1 + q) 7 / 2 
(1 + x)(l — 3x) 

X x 2 (l-x ) 3 ’ 

{2 + (1 - x)q} 2+{ - 1 - x '> q 


(16) 


[{1 - (1 + 2a;)g}/2]{ 1 -( 1 + 2 *M/ 2 {3(l + q)/2 } 3 d+9)/2 
x{x"(l-^) (1 ^ ) }- 2 «. 
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(17) 











The upper value of the range of integration is given by 

1 for 0 < q < | 

^ ^ \<q< I- 

We shall first make the Gaussian approximation for ^t[x, q) N as a function 
of x, and obtain 


fN q ^N~i f dxcj> T {x, q)<S> T {x, o, q) N e -" a °^ x ~^ 2 . 

Jo 

Here, the peak position Xq is determined by solving the equation 

ds T (x, q) _ 
dx 

where we define the asymptotic entropy density as ST{x,q) = — logf'h-r^, q)}, 
and we obtain xq = (1 — q )/{3 — q). We also write 


ao{q) 


d 2 s T {x,q) 

Ox 2 la:=:co, 


which turns out to be (3 — q ) 3 (2 — q)q /(1 — q) 2 /6 . 

Since xq is smaller than x max for any q within the range [0,1), the integration 
over x is dominated by the contribution around the peak xq. After carrying out 
the integration, we repeat the Gaussian approximation for &T{xo(q),q) N as a 
function of q , and obtain 


/at, q 


N 2 q 2 


(2 - go ) 1/2 


(1 + < 7 o ) 7 / 2 (3 — go ) 1 / 2 


$ r{x 0 ,qo ) 


N 


ex P{-y^o {q - go) 2 }, 


where qo is determined by the eciuation 


ds T (x 0 (q),q) _ 
dq 

which we solve numerically obtaining q$ ~ 0.2067, and where we write 

A _ d 2 s T (x 0 (q),q)! 

— Qq 2 19 = 90 ’ 

which is about 2.682. Finally, the asymptotic distribution appears precisely in 
the same form as other cases: 

fN, m -*rn^N~i exp(-^m 2 + fim - A N), (18) 

with fj, = A 0 <7 o, and A = A 0 q%/2 - log[$ T (a:o(go), go)]- 

It is worth noticing that despite the apparent difference of the Tutte distri¬ 
bution function from others, all the distribution functions in asymptotic regions 
fall in the same universal form with the same critical exponents m 2 N~ 2 . This 
suggests that even a simple product of generalized binomial coefficients, which 
are supposed to reflect the hierarchical structure of the manifold, can produce 
the universal form of the distribution function in the asymptotic limit. We shall 
rely on this observation for obtaining the 3-dimensional distribution function in 
the next section. 
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3 A mathematical model of the distribution 
function for the 3-disk 

So far no recursion equation of the generating function is known for three di¬ 
mensions, and we would like to propose a mathematical model inspired by the 
knowledge we have gathered in 2-dimensional discrete manifolds. We shall con¬ 
struct the 3-disk with a simple S 2 boundary from the 3-sphere. 

3.1 construction of the discrete 3-disk 

The 3-dimensional sphere is prepared as a simplicial complex with the total 
numbers of i-simplices being A/j (i = 0,1, 2, 3). The Euler-Poincare relationship 
for the S 3 topology gives 


Mo ~ Afi + M 2 — M 3 — 0, (19) 

together with the dimensional constraint, 

A f 2 = 2M 3 . (20) 

The 3-dimensional disk may be obtained from the 3-sphere by choosing one 
vertex among Mo, and removing all the tetrahedra sharing the vertex. We call 
the removed complex as the cap, and the remaining disk as the target manifold. 
They both have the D 3 topology with a common S 2 boundary. The number 
of tetrahedra x of the cap is called the coordination number of the vertex. 
Numbers of i-simplices for the cap and the target manifold are written as {C;} 
and {Ni}, respectively. For the target disk TV, are decomposed into sums of 
those belonging to the boundary TV,; and its complement TVj satisfying 

Ni = Ni + Ni (* = 0 , 1 , 2 ), 

which are often called the external and the internal i-simplex, respectively. Since 
the original 3-sphere is divided by the common S 2 boundary into two 3-disks, 
we have the relationships 

Mi = Ni + Ci-N, (* = 0 , 1 , 2 ), 

M 3 = N 3 + C 3 . (21) 

According to the construction procedure the cap variables are given to be 


C 0 = N 0 + 1, Cr = TVo + TV 1; C 2 = TV, + TV 2 , C 3 = 1 V 2 . 
By adding the two topological conditions of the S 2 boundary: 


( 22 ) 


TV 0 - TVi + TV 2 = 2, 2TV! = 3TV 2 , 


(23) 


we count altogether 15 independent relationships among 18 variables. 





Choosing No,Ni, and N 2 as the independent variables, other variables rep¬ 
resenting the target 3-disk are expressed as 


N 0 

= N 2 - 2N X + 3 N 0 + 4, 


N, 

= 3 N 2 - 5W + 6 No + 6, 


n 2 

= 3 N 2 - 4 Ni + 4 N 0 + 4, 


n 3 

= N 2 - JVi + N 0 + 1, 

(24) 


and 


N 0 

= N 2 - 2Ai + 2 No + 4 , 


Ni 

= 3 N 2 - 6 N 1 + 6 N 0 + 6, 


n 2 

= 2 N 2 - 41Vi + 4 N 0 + 4. 

(25) 


Let us consider the practical construction procedure of a possible 3-disk 
(No, Ni, N 2 ). We start with an elementary tetrahedron (0,0,0), and either attach 
or detach one tetrahedron to or from the target disk one by one just like a crystal 
growth through the diffusion limited aggregation. There are 3 ways to attach 
a tetrahedron corresponding to the number of neighboring faces to be glued as 
shown in Figures 1, 2, and 3. Each Figure shows a pair of circles. The left 
circle and the right circle represent the S 2 boundaries before and after gluing 
a tetrahedron, respectively. Solid lines and dotted lines represent the external 
and the internal edges, respectively, while the dashed lines are also the external 
edges, which are hidden behind the attached tetrahedron. These gluing and 
removing procedures are called moves, which we represent as A 4i(i = 1,2,3), 
and their inverses as A 4i- They are sometimes called up-moves and down-moves, 
respectively. Application of these three up-moves on a target 3-disk changes of 
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Figure 2: A^-move 



Figure 3: A^-move 
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variables A r ,. and Ni, as listed in the following Table, 


moves 

A N 2 

A No 

AN 0 

A N 3 

A N 2 

A No 

AN 0 

an 2 

Mo 

1 

0 

0 

1 

3 

3 

1 

2 

m 2 

2 

1 

0 

1 

2 

1 

0 

0 

CO 

3 

3 

3 

1 

1 

1 

0 

0 

-2 


where we add AIV 2 in the list for the later use. 

Let us first show that any possible 3-disk complex can be constructed by 
appropriate combinations of these three types of up and down moves. Writing 
a net number of the A4;-move as mj, i.e. number of A4i-moves minus number of 
A4.;-moves, the target 3-disk acquires the following numbers of simplices: 


No 

No 

N 2 

N3 


and 


also 


No 

No 

N 2 


No 

No 

N 2 


4 + mo, 

6 + 3toi + m 2 , 

4 + 3toi + 2 m 2 + m 3 , 

1 + rrio +m 2 + m 3 , 


= 4 + TOI — to 3 , 

= 6 + 3toi — 3m 3 , 
= 4 + 2 toi — 2 m 3 , 


= rn 3 , 

= m 2 +3m 3 , 

= mo + 2m 2 + 3to 3 . 


(26) 


(27) 


(28) 


Inserting mi(i = 1,2,3) written in terms of Nq,No,N 2 obtained from the last 
three equations: 

mo = N 2 — 2No + 3 N 0 , m 2 = No - 3 N 0 , m 3 = N 0 , 


the remaining 7 equations give the relationships, eq. (24) and (25). This implies 


that the three types moves are sufficient to construct a simplicial complex with 
any possible number of variables, {Ni, Ni, Ni}. However, this does not neces¬ 
sarily mean any possible complex can be constructed by some combinations of 
those 3 types of moves. 

Secondly, let us show that it is indeed possible by searching a path starting 
from a certain possible complex, and choosing appropriate down moves to arrive 
at an elementary tetrahedron. Then, the inverse path gives one of the actual 
construction of the desired complex. We start with a D 3 complex (No, No,N 2 ), 
which has N 2 external faces. Each external face is sorted to the following 3 
types: a) one of 3 faces of a 3-simplex attached to the main body by sharing 
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one face such as the complex created by an Adi-move. b) one of 2 faces of a 
3-simplex glued to the main body sharing 2 faces such as that created by an 
Af 2 -move. c) the face of a 3-simplex glued to the main body sharing 3 faces such 
as that created by an Af 3 -move. It is possible any time to remove the 3-simplex 
whose external face belong to the a-type. It is not true when the face is either 
b- or c-type, in case when the down-move results non-allowed triangulation such 
as creating a pair of 3-simplices which are attached through a vertex or an edge, 
or creating a vertex or an edge whose neighboring simplices do not form the D 3 , 
or the semi-Z ? 3 when it is on the boundary. 

We begin with detaching all those simplices whose faces belong to the a- 
type. Then the remaining external faces belong either to the b- or the c-type. 
We select any one of simplices which are allowed to remove. After the down- 
move we again check all the external faces and remove all those simplices whose 
faces belong to the a-type. Repeating these procedures mi + m 2 + m 3 times we 
eventually arrive to an elementary tetrahedron. Inverting these moves starting 
from a tetrahedron we can create the desired complex. 

From the above construction procedures we recognize that there exist three 
hierarchical structures created by combinations of the 3 types of moves: i) 
applications of the Afi-moves create complexes without any internal vertex and 
internal edge. Then, all the internal triangles created by this type of moves 
have their associated three edges and vertices on the boundary like a branched 
tree, and thus these complexes belong to one-dimensional sub-manifolds, ii) 
Combinations of the Adi- and Ad 2 -nroves create complexes having no internal 
vertices. In such a complex two end-vertices of all the internal edges must lie 
on the boundary surface, that constrain the thickness of the 3-disks being the 
minimum. It characterizes the manifold to be 2-dimensional like, iii) Applying 
all the three types of moves, the resulting 3-disks are the 3-dimensional generic 
manifolds. 


3.2 the model distribution function 

We assume maps are rooted. A rooted map of the 3-disk is obtained by speci¬ 
fying an external triangle as the rooted triangle, and selecting one of the three 
vertices and an arrowed edge pointing to the vertex under the right-hand screw 
rule directing outward as the rooted vertex and the rooted edge, respectively. 
It is sometimes more convenient an d ph ysical to choose N 2 as the independent 
variable rather than N 2 through eq.(25), which we shall employ hereafter. 

The distribution function of the rooted 3-disk is expressed by m , where 
we define 

n = N 0 , l = N u m= |(JV 2 — 4). (29) 

By removing a cap with the coordination number x from the S 3 manifold , we 
obtain the target disk with 


Ah 


= ^- 2 - 3 ’ 
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together with 


( 30 ) 


Ni = U 0 +U 3 -2x-2, 
N 2 = x, 


N 3 =AT 3 -x. (31) 

Therefore, once we obtain the coordination number distribution of the S 3 com¬ 
plex, P^f 0 ^ 3 (x), we can calculate the distribution function for the target 3-disk, 

fnlm, aS 

fnj,m = PMo,M 3 {x)- (32) 

Corresponding to the three hierarchical structures there are three types of 
distribution functions. When we apply only the Ali-move on the manifold 
starting with a single tetrahedron, we can produce the restricted distribution 
function m corresponding to the 1-dimensional sub-manifold. Maps created 
by the Ali-moves are expected to be identical to the 2-dimensional quadrangu¬ 
lar maps with zero internal vertex, which have the branched tree distribution 
function <?o,m of eq. ([7]), constructed by tetravalent vertices: 

f D 3 _ o( 3 m + 2 \ 

/ 0,0,m 2m + 3/ 

Here, we restrict ourselves considering only spacial degrees of freedom without 
any matter field coupled (pure gravity ) for simplicity. 

By including the XI 2 -moves we shall obtain maps containing non-zero N 3 
and N 2 . They give the distribution function of the form rn . We expect 
from observations in the previous section that when the maps contain the 1- 
dimensional sub-manifold, the distribution function will be given by a product 
of two generalized binomial coefficients similar to gi, m eq. (J6]) . By making use of 
all types of moves we expect to construct a generic 3-dimensional manifold. 

Now we give a naive but natural ansatz for the distribution function of the 
3-disk, which possesses the double hierarchical structure, by a product of three 
generalized binomial coefficients, 




f 1 

J n,l,m 


= XA n 


x C n 


a\m + a 2 \ z / b\m + b 2 l + b 3 
to, a 3 m + a^J \l, bim + b^l + bg 
cirri + c 2 l + c 3 n + C4 
n, C5to + cqI + crn + c$ 


(33) 


In this expression we count total 22 parameters, A, B , C, X , {cq}, {&;}, and {c^}, 
among which we are going to eliminate by imposing various conditions required 
either from theoretical or numerical reasons. From the experience in the 2- 
dimensional exact distribution functions, the parameter sets at,bi, Ci is expected 
to take either positive integer or half-integer values. We consider here limiting 
within integer values for simplicity. 

We should remind that the model distribution does not respect details of 
practical construction rules, such as the n-conectedness, or acceptance of de¬ 
generate graphs. Relying on the observations in the 2-dimensional exact dis¬ 
tribution functions, we expect it is only in the asymptotic region the model 
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distribution function reveals the universal form, which is free from constraints 
put in practical construction procedures. 

One of the most basic theoretical requirement, which is not proven mathe¬ 
matically but seems to be true numerically, is the exponential boundedness of 
the number of independent geometries by the total number of 3-simplices N 3 . 
According to the asymptotic limit of the generalized binomial coefficients, the 
most diverging factor N n is required to disappear. It is satisfied by imposing 
three conditions, 


ai - 1 - a 3 + bi — 64 + Ci - c 5 = 0, 

62 ^ 1 — 65 + C2 — C6 = 0, 

c 3 - 1 - c 7 = 0. (34) 


Hierarchical structure of 2-dinrensional sub-manifold requires that the dis¬ 
tribution function of the 2 -dinrensional sub-manifold 


/ aim + a 2 \gif &i m + M ffe \ 
\m, a 3 m + a^J \l, b^m + b 5 l + b e J 


{ C\m + c 2 l + C 4 
\ 0 , c 5 m + c 6 l + c 8 


should be written as a product of two generalized binomial coefficients. It is 
achieved by asking 

b\ = c 5 , b 2 = c 6 , 63 = c 8 , (35) 


resulting 


/ ai m + a 2 \ ; / cim + c 2 l + C 4 \ 

\m, a 3 m + < 3 , 4 / \l, bim + 65 / + bej 


When no matter degree of freedom is involved, we ask the distribution func¬ 
tion of 2-dimensional sub-manifold has the same critical exponent as the 2-disk 
distribution function: 

C 4 — be = —2. (36) 

Similarly, the distribution function for the one dimensional sub-manifold, 


fD _ y An 

I 0 , 0 ,m — 


a\m + a 2 
m, a 3 m + 04 


/ cim + c 4 \ 
\ 0 , 64771 + be)’ 


is required to be that for branched trees constructed by tetravalent vertices 
given by g 0 ,m eq. §■ This gives 

ai = 64 , a 2 = b e , (37) 


and 

X = 3, A = 1, ci = 3, C 4 = 2, a 3 = 2, 04 = 3. (38) 

At the end we obtain the 3-disk distribution function for pure gravity as 


f D 

J n.l.m 


X 


/ oiTO + 4 \ hm + b 2 l + 5 \ 
\rn, 2m+ 3J \l , aim + 65 / + 4/ 
m / 3m + (65 + 1)/+ c 3 n + 2 \ 
\n, bim + b 2 l + (c 3 - l)n + 5/ ’ 


(39) 
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where the parameter 63 is set to be 5, so that C 4 — Cs = 2 — 63 being — (manifold 
dimension ), which has been the case in 1 and 2 dimensions. In any way 63 will 
not appear explicitly in the asymptotic limit. We still count 7 free parameters at 
this moment left undetermined, which may be fixed by comparing to numerical 
simulations in the following sections. 


3.3 asymptotic distribution 

In order to study the universal behavior of the distribution function which are 
expected to appear in asymptotic regions, we define new variables similar to the 
2 -dimensional cases corresponding to the independent variables n,l,m: 

p = n/N, r = l/N , q = m/N, (40) 


where 


N = N :i - 1. 


From eq.(24) and (251 there is a relationship, 


q + r — p = 1 . 


(41) 

(42) 


We shall select (N, p , q) as the new set of independent variables when we discuss 
the asymptotic limit. 

By making use of the asymptotic formula of the generalized binomial coeffi¬ 
cient eq. the model distribution function of the 3-disk in asymptotic regions, 
which we express as ( p,q ), is written within an overall constant by 


SnM ~ N-*l 2 qV 2 <l>{p, q )*b, q ) N , 


(43) 


where 


(pi q + b 2 r) 11 / 2 (3q+ (fr 5 + l)r + c 3 p ) 5/2 

9 r 1 / 2 (a \q + b 3 r) 9 / 2 p^-^ibiq + b 2 r + (C 3 — 1 )p) 11 / 2 ’ 


(44) 


and where 


Hp,q) = (ar/4) q B r Cr 


(biq + b 2 r) b ^+ b * r 


r r (aiq + 6 5 r) ai9+f>5r ' 

{3g + (65 + 1 )r + c 3 p} 3<!+{b,i+1)r+C3P 
p p {biq + b 2 r + (c 3 - l)p} b i<?+ b 2 r+(c 3 -i)p' 


(45) 


4 Phases of the stationary distributions 
and comparison to numerical results 

We now wish to check our proposition for the distribution function of the dis¬ 
crete 3-disk by comparing to numerical results. Until now there are number 
of numerical studies available US [H ES EH US] for pure gravity with S 3 
topology, and one result of D 3 [27] . We shall look for appropriate values of 
parameters left undetermined in the last section from the numerical data. 
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4.1 S 3 topology 

We can study statistical properties of the discrete 5' 3 -manifold through the 
canonical partition function 

= E eVNo ( 46 ) 

T(S 3 ) 

= E/C^, 

A/o 


where sums T(S 3 ) are taken over all the independent triangulations with S 3 
topology. The distribution function /jy - 0 can be obtained from the D 3 dis¬ 
tribution function with the least boundary area, which is created from the S 3 
simplicial complex by removing one of A 3 3-simplices, as 

/Jv'oA3 = fn,l,o/^i ( 4 ^) 

where 

n = A/o — 4 , l = A/o + A /3 — 6 . 

In the asymptotic limit we get 

/Az-oW 3 /n(p) ~ N-^(p)(BD) n (BCDY n Hp) n , (48) 


where N = A /3 — 2,p = n//V, and D = b^ 2 /b$ . After separating those terms to 
be included in either the normalization constant or the cosmological constants 
in <j>{p, 0 ) and 4>(p, 0 ), we have 


4>(p) 

Hp) 


1 

pi?*’ 

r (b 2 -b 5 - l)r {(^ + + C3p }(65 + l)r+c 3 p 

pP {/) 2 r + (c 3 - l)p} b ^+(c 3 -i)p ’ 


(49) 


with r = 1 +p. We note the analytic part in <f>(p, 0) is set to a constant because 
of $>(p) N factor which dominates the sum over A/o in the large N limit. 

Asymptotically the partition function can be written in the integral form as 



N [ PmaX dp f s N \p) e vpN . 

Jo 


(50) 


Upper value of the p-integration p max 
N —> 00 : 


is determined from eq. (31) in the limit 


P 

q 


V-Z/2 
i-£ ’ 
£ 

2 ( 1 - 0 ’ 


(51) 
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where £ = x/Mj, and 77 = A/0/A/3. Since the physical range of q is [0,1), £ 
varies within the interval [0,2/3), and p takes the value in ( 3r) — 1 , 77 ] . When 
q ~ 1 (i.e. £ ~ 2/3), most of the tetrahedrons have two neighbors to form a 
linear chain with its two faces attaching to neighboring tetrahedrons, while the 
remaining two faces left as external faces. This branched tree structure does not 
contain internal vertices, and the value p = 0 at this limit fixes the maximum 
value of 77 to be 1/3, as conjectured by Walkup [55] . The upper limit of p is 
fixed at the maximum value of 77: Pmax = 1/3. 

One of the important findings of DT simulations is the nature of phase 
transition, which turned out to be of first order. Property of the phase transition 
can be studied through the free energy density defined by 

s S n(p) = ^l°g[/jv(p)]- ( 52 ) 

Asymptotically it converges to 

S N (P) =*■ 


The second derivative of the free energy density, 

( 2 ) / , ^ c 3 (c 3 - 1)(1 + b 5 - c 3 )p + b 2 { 1 + b 5 )c 2 r 

SN P p(p + 1)((1 + b§)r + C 3 p)(b 2 r + (c 3 — l)p) ’ 


d 4. fr. ji+fcs 

lo s[-^-] (p-> 0), 

(1 + 65 + c 3 ) 1+b5+C3 

lo e [ (i2 _ 1 + c 3 ) h-i + =,]J>(i> 


(53) 


behaves as —1/p for small p, and as p gets larger it increases and cross zero once 
at p = pi with 


Pi 


1 

( 65+T — bj) c 3( c 3 — 1) — 1 


(55) 


When pi lies below 1/3, the distribution function has two peaks, one at 1/3 
and another between 0 and p\, for an appropriate parameter range of v. In this 
case the system makes 1 st order phase transition at the critical v, where heights 
of the two peaks are equal. Numerical simulation shows that the average Nq 
has a jump at v = v c ~ 4.0, with a small gap of <5 < p >~ 0.03 [24j, and we 
shall choose the integer parameter triplet (&2 7 ^ 5 ; C 3 ) in the region where p\ is 
less than but close to 1/3. We note that the system shows second order phase 
transition at pi = 1/3, and it becomes cross over transition as p\ gets bigger 
than 1/3. 

In order to fix parameters from the DT simulation results we need to make 
direct numerical integrations for the canonical partition function and other ex¬ 
pectation values. Among the unknown parameters there is no observable to 
determine the values directly except B 1 and C , which are related to the cosmo¬ 
logical constants. We shall fix the parameter triplet (& 2 , & 5 , c 3 ) which control the 
nature of phase transition through p\, by asking 1) they are positive integers 
and choose as small as possible just for simplicity, 2 ) they give p± reasonably 
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Figure 4: < p > (v). N = 100. v c is fixed at 3.6. 


close to 1/3 in order to make the jump in the average Nq small, and 3) they 
give correct susceptibility peak pattern in the 2 -dimensional parameter space 
of the D 3 topology, which will be discussed in the next sub-section. At first 
we fix b 5 to be 3, after checking that b 5 = 1, and 2 do not satisfy the above 
condition 3. With the fixed b§ values of 62 and C 3 are bounded below by the 
requirement p\ < 1/3. Next, we choose 62 = 6 , which is the smallest possible 
value, and lastly C 3 to be 8 in order to make pi close to 1/3. With these pa¬ 
rameters the expectation value < p > (~ A/0/A/3) and the finite size scaling of 
the A/"o-susceptibility are shown in Figures 4 and 5, respectively, which clearly 
exhibits the characteristics of first order phase transition. 

Two constants B and C are fixed from the position u c of the jump in the 
average value < Nq >, and the asymptotic entropy density po obtained from 


the canonical partition function eq. (50), which we parametrize as 


^log(Zf ! [0])~^ 0 -(3-7) 


log N 
N ’ 


(56) 


at a point v = 0. In Figure 6 we show log(Z^ 3 [0])/iV with a fixed po ~ 2.0, 
which gives 7 ~ 0.35. This parametrization differs from numerical simulations. 
There were two numerical results available for the parameter po'- i n Ref.|23j it 
was 

^l°g(zf i [0])~Aio-/3 i ^, 

with po ~ 2.06, a ~ 0.32 and [3 ~ 3.90 , while in Ref.[22] it was 

V log(z »[o])~»- /J toiw 

with po ~ 1.62 and f3 ~ 1.44. Our parametrization intends to give a meaning 
to the parameter 7 to be the string susceptibility of the discrete S 3 manifold. 
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Figure 5: Finite size scaling of x( v n) for N = 100,200,300,400,500. vn is 
chosen at the peak position of x(iz).The line is for guiding eyes. 



Figure 6 : Numerical fit of log {Zf^[v])/N by n o — (3 — 7 ) lo fj N , with /x 0 = 2.0 
and 7 = 0.35 at u = 0. 
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Numerically po and 7 may depend on v. In Figures 7 and 8 we show the v 
dependence of p 0 and 7 obtained by numerical fitting of log (Zfj [v])/N for N 
between 1000 and 10000, respectively. They are almost constant up until the 
transition point v ~ 4, and after that they behave differently. This is because 
of the limit in the integration over p at p max = 1/3. When v > r/ c , the peak at 
p = 1/3 dominates in the integration. As for the constant i/ c , there exist also 
two simulation results: ~ 4.0 determined from the S 3 DT-simula tions [231 2U, 
and ~ 3.6 observed in the susceptibility map of the D 3 DT-simulation [ 25] . 
These facts pretend us existence of the strong finite size effect and probably the 
systematic errors depending on numerical methods. 

4.2 D 3 topology 

Since the current study is intended for a construction of the 3-dinrensional Uni¬ 
verse, which stands alone in the non-framed ’space’, there is no meaning to 
specify a root in the map, and the distribution should be counted for non-rooted 
maps. Transformation of distribution functions from rooted to non-rooted D 3 - 
maps is approximately performed through dividing the Z3 3 -distribution func¬ 
tion by the number of possible root sinrplices, which equals to the number of 
boundary 2-sinrplex N?- If a map has spacial symmetries, such as rotational, 
reflectional, etc. this simple rule cannot be applied. However, such symmetric 
maps are rare especially when N is large, and we expect they will not give much 
effects on the distribution function in the asymptotic limit. We shall write the 
approximate distribution function of non-rooted maps as 

=—f D3 

Jn.l.m _ J n.l.m * w ' / 

’ ’ 171 ’ ’ 

Statistical properties of the discrete D 3 -nranifold are studied through the 
canonical partition function 

Z% 3 [»A= E e- SM , (58) 

T(D 3 ) 

where sums T(D 3 ) are taken over all the non-rooted maps with the D 3 topology. 
Here the action for the D 3 topology is given by 

S[»,v}=pN 2 -vN 0 , (59) 

In the asymptotic limit the partition function is written in the integral form as 

, r 1/3 rq ma *{p) .... 

Z° \nM =N 2 dp dq fg (60) 

Jo Jo 

where p = n/N and q = m/N as before, and where 

f£ 3 {P,<l) = fnlm- 
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The distribution function f® ( p , q) is defined from the original distribution func¬ 
tion by replacing l to N using the topological relationship N = l — n + m. 
For a fixed p, eq.((5l]> gives q ma x{p) = 1-3 p. 

The model distribution function for non-rooted D 3 - maps in the large N limit 
is written as, 


InM ~ N- 7 / 2 q- 1 / 2 $(p,q)B N (BCr N (a a 1 l /A/By N ^(p,q) 


N 


(61) 


where 

= rl /2pl/2 » ( 62 ) 

and 

q (a 1 -3)q ( biq + b2r jb iq+ b 2 r + ^ + + C3? ,)3g+(b 5 + l )r+c 3 p 

yP’Q)- r r pP ( aig + b 5 r) b ^+ b ^ r {biq + b 2 r + (c 3 - l)p) & i9+f>2r-+(c 3 -i)p' 

(63) 

In the above equations r should be replaced by 1 + p — q. Analytic part of 
<t>(p, q) is set to a constant in the same reason as the S 3 case, and also in order 
to suppress the N dependence in numerical integrations with finite N. We need 
to make the double integrations for the patition function and other expectation 
values, which is time consuming. For getting a better convergence and to save 
computer time we set N as small as possible, and choose N = 30 in the following 
calculations for the susceptibility landscapes in the two parameter space. 

To carry out the numerical integrations, we need to fix additional two in¬ 
tegers for <22 and b± for <1 ){p,q). Starting with a\ = 1 and &i = 1, we find ai 
should be larger than 3 under the similar stuation as the condition on b$. After 
fixing a± = 3 we search appropriate b\ starting from 1. Until b\ reaching to 
5 the susceptibility peak range lying vertically is tilted, and we finally select 
(ai, &i) = (3, 5). We need to adjust the parameters B and C in order to fit the 
horizontal range of susceptibility peaks lying along v ~ 3.6, instead of the S 3 - 
sumulation value of 4.0, and also p c = 1.95, which lyes between two numerical 
results 1.6 and 2.0. The agreement between the numerical result |25j and our 
model turned out to be unexpectedly good as we see in Figure 9 for the average 
N 0 , and Figure 10 for the 7V 0 -susceptibility. 

We also calculate the average values of p ~ N 0 /N, and q ~ N2/N as shown 
in Figures 11 and 12, respectively. They show clearly the existence of three 
phases reflecting the hierarchical structure of the U 3 -manifold: i) The platou 
in < p >~ l/3(blight area) correspond to the 3-dimensional generic manifold, 
ii) The blight area in < q >~ 1 represents the ID sub-manifold (branched 
tree) phase, and iii) the remaining area, which is the dark areas in both < p > 
and < q > is considered to be the 2-dimensional sub-manifold (< r >~ 1 ). 
These three areas are surrounded by susceptibility peaks determined by the 
DT-simulations. 
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Figure 9: < N 0 > (/x, v)/N. N = 30. 



Figure 10: 7V 0 -susceptibility. Red dots are taken from |2S] with dashhed lines 
for guiding eyes. 
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Figure 11: < p > {n,v). Hights of the blight platau is ~ 1/3 and the dark 
bottom ~ 0. 



Figure 12: < q > (n, v). Hights of the blight platau is ~ 1 and the dark bottom 
~ 0 . 
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5 Discussions and Summary 

For predicting the distribution function of the discrete 3-disk we have made use 
of the common structural property in the exact forms of certain 2-dimensional 
disks. We rely on the following two observations. Firstly, when a discrete 2- 
disk manifold contains 1-dimensional sub-manifold the distribution function is 
written in a simple product of two generalized binomial coefficients. Secondly, 
even if a manifold does not belong to above mentioned class, such as the Tutte 
distribution function, the distribution function will show asymptotically the 
same universal form with the same critical exponent as the product type. In 
the 3-dimensional case the true distribution function of a manifold constructed 
by the DT may not be written in a simple product form. However, we know 
that the discrete 3-disk constructed by the DT method possesses hierarchical 
substructure which persists in the asymptotic regions, and we expect a simple 
product distribution function of generalized binomial coefficients will show the 
same universal asymptotic behavior as that obtained by the DT. 

In oder to check above prediction we compare our formula of the discrete 
3-disk distribution function with numerical results obtained by the DT simula¬ 
tions. We choose a set of parameters which fulfill exponential boundedness of 
the state volume, and show the hierarchical structure of sub-manifolds and 1st 
order phase transition, which is one of the significant feature of DT simulation 
in S 3 manifold. In general, our predicted distribution function shows the pos¬ 
sibility of three types of transitions, either crossover, 1st order, or 2nd order, 
depending on the choice of the parameter set. Since 2nd order phase transition 
appears at a point on a line in the parameter space where phase changes from 
crossover to 1st order, we understand 1st order transition has more probability 
of appearance. There will be a possibility to obtain 2nd order phase transition 
if we put appropriate matter fields on 3-dimensional discrete space. 

Phase structure of the 3-dimensional manifolds is another important point 
to check. It can be visualized by the expectation values of p and q in the two 
dimensional parameter space (/x, v). We observe three characteristic regions: i) 
for large v with small /x, a plateau with the hight < p >~ 1/3, which corresponds 
to generic 3-dimensional manifold, ii) for large /x, a plateau with < q >~ 1 which 
exhibits the 1-dimensional sub-manifold property, and iii) small v and /x where 
< p >, < q 0, which is expected for the 2-dimensional sub-manifold. Three 
phases are well-separated by two walls in the susceptibilities: one in \ p which 
separate the generic 3-dimensional disk and the 2-dimensional sub-manifold, 
and another in \q which divides the 1-dimensional sub manifold from others. 
Over all view of phase structure is quite similar to the numerical result [25] . 

Our primary task of achieving an analytical form for the distribution function 
of 3-disk is to apply it for the statistical evolution equation of the 3-dinrensional 
Universe. The ’free energy’ log]/] acts as the force field in the master equation 
for the probability of the 3-dimensional Universe with the volume V = the 
surface area A = fly/V/, and the ’dilatation’ D = dgN 0 along the ’Markov-time’ 
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t. When we define the physical time t by 


t = 



1 d,V{r) 
A(t) dr 


dr, 


(64) 


so that it satisfy the rate of change of 3-volume to be dV{t) = A(t)dt , the 
relationship between the physical time and the Markov time is given by 


/ tV 1 - ^ ^ i 

1 log T if / = 1 


for a solution of the master equation with the asymptotic behavior, V ~ r and 
A ~ t k . If we ask the physical time to posses the physical dimension 1, then the 
Markov time r will have a dimension 1/(1 — k), and thus the body (V) and the 
face (A) carry the dimension 1/(1 — k) and k/(1 — k), for the case k/ 1, respec¬ 
tively, while for the case k = 1 the space dimensions are infinity. We recall in 
the 2-dimensional Universe the statistical evolution equation had three asymp¬ 
totic solutions, k = (1,1/2,0), and thus the face dimension varies (oo, 1, 0). We 
expect for the 3-dimensional Universe there will be a k = 2/3 solution in the 
master equation, so that the dimensions V and A are 3 and 2, respectively. We 
are now underway to obtain solutions for three dimensional evolution equation. 
We also focus our investigation on the effect of phase transition in the evolution 
of the Universe. Extension to the 4-dimensional case is straight forward, and it 
will soon appear in separate articles. 

We would like to add some comments on the colored tensor model which has 
been claimed to be a generalization of the matrix model for higher dimensions. 
It is then regarded to be an analytic realization of the dynamical triangulation 
of space-dimensions higher than two. In the model the dominant contributions 
come from graphs constructed with an elementary block called melon, which is 
formed by a pair of 3-simplices glued along 3 faces. The generating function of 
rooted melon graphs is known to be 

T( A) = 1 + A 2 T(A) 4 


where A is the vertex coupling constant in the dual graph of the 3-dimensional 
triangulation. This equation can be interpreted as the generating function of 
branching trees made by 5-nary vertices, which is known to give the distribution 
function, 

4m A 

m, 3 m + 1/ 

where m = (N 2 — 4)/2 = N. Asymptotically the number of graphs is Tjv ~ 
N -3/2 e » 0 N _ ( 1 / 2 ) log(256/27), which is about 1.12. On the other hand 

in the DT simulations it is either 1.6 or 2.0, which considerably surpasses the 
nrelonic graphs. We also remark that in the standard DT simulations the melon 
graph is excluded as degenerate graphs. From the distribution functions T m we 
can observe that the nrelonic graphs belong to the one dimensional universality 
class. Recently, a new type of tensor model is introduced(30] with the critical 
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property showing the 2-dinrensional character in the double scaling limit. We 
expect that another new tensor model which exhibit a genuine 3-dinrensional 
phase transition will appear, so that the transition properties can be compared 
with our phenomenological model quantitatively. 
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